function [x,path]=prox_nag(prox,init,theta)
iter=1e5;

v1=init;
dim=length(init);
v2=prox(v1,1);

path=zeros(dim,iter);
for i=1:iter   
    % a=1;
    a=1/(i+3)^(8/9)/4;
    x=(1+theta)*v2-v1*theta;
    v1=v2;
    v2=prox(x,a);
    path(:,i)=x;
end
